Nonlinear Rashba Model and Spin Relaxation in Quantum Wells 
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We find that the Rashba spin sphtting is intrinsically a nonlinear function of the momentum, 
and the linear Rasha model may overestimate it significantly, especially in narrow-gap materials. A 
nonlinear Rashba model is proposed, which is in good agreement with the numerical results from the 
eight-band k-p theory. Using this model, we find pronounced suppression of the D'yakonov-Perel' 
spin relaxation rate at large electron densities, and a non-monotonic dependence of the resonance 
peak position of electron spin lifetime on the electron density in [lll]-oriented quantum wells, both 
in qualitative disagreement with the predictions of the linear Rashba model. 

PACS numbers: 71.70.Ej, 72.25.Rb, 73.21. Fg 
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Recently there has been growing interest in the field 
of SpmtromcsP, 12, 13, which explores the electron spin, 
in addition to the electron charge, to realize new func- 
tionalities in future electronic devices. 0, IE ^ promis- 
ing approach implementing such spintronic devices is to 
utilize the Rashba spin-orbit interaction caused by struc- 
ture inversion asymmetry in quantum wells (QW's),|3| 
which can be controlled by gate voltagesj^ as well as 
band structure engineering. JdJ Approximate analyt- 
ical expressions based on second- or third-order pertur- 
bation theories [Tol Ull IT^ suggest that the Rashba spin 
splitting (RSS) is a linear function of the in-plane wave 
vector fc|| . This linear Rashba model has been widely 
used to investigate the various spin-related properties of 
low-dimensional semiconductor structures, e.g., electron 
spin relaxation^E llH ll^ ll^ an d the ne wly discov- 
ered spin Hall effecijlilSnifl Ill^J^ . How- 
ever, recent numerical calculations |l2l l25l bsj show that 
the RSS in certain semiconductor QW's deviates from 
the linear behavior at large fc|| , although the underlying 
physics remains unclear. Since the linear Rashba model 
is still widely used by the mainstream researchers, it is 
important to explore the underlying physics beneath such 
deviation and, if necessary, check the validity of the linear 
Rashba model. 

In this work, we find from analytical derivation from 
the eight-band k • p theory that the RSS is intrinsically 
a nonlinear function of the wave vector, which is caused 
by the weakening of the interband coupling with increas- 
ing kinetic energy of the electron in the conduction band. 
We show from numerical comparisons that the deviation 
of the linear Rashba model could be surprisingly large 
at large wave vectors, especially for narrow-gap materi- 
als. These facts substantiate the necessity for a nonlin- 
ear Rashba model. We propose such a model, which is 
in good agreement with the numerical results from the 
eight-band k • p theory for various QW's. This nonlin- 
ear Rashba model would lead to a series of significant 
modifications to the various spin-related properties of 



the electron, most of which has been investigated based 
on the linear Rashba model. For example, we find pro- 
nounced suppression of the D'yakonov-Perel' (DP) spin 
relaxation rate (SRR) at large electron density, in qual- 
itative disagreement with the prediction of the linear 
Rashba model. 15,, IfiJ The resonant enhancement of the 
electron spin lifetime in [lll]-oriented quantum wells|l7j 
also exhibits qualitatively different behavior from the lin- 
ear Rashba model. The values of the SRR obtained by 
the two models differ by up to several orders of magni- 
tude. Further surprising consequences are expected when 
this nonlinear Rashba model is applied to other fields. 

For [001]-oriented heterostructures, the eight-band 
Hamiltonian[2^ H = Hk + V , where V = eFz is the ex- 
ternal electric field induced potential and Hk is the eight- 
band envelope function Hamiltonian. Neglecting the off- 
diagonal elements in the valence bands and eliminating 
the valence band components of the envelope function, 
the effective conduction band Hamiltonian is obtained as 



H, 



e//(fc||) = E^{z) + V{z) + k— k + ao(z)(fc|| x e,) ■ a, 

(1) 

where ao = / {&ma)d^{z) / dz, 7(2) = Ep[l/Uih{z) - 
\/Uso{z)]. Ep = 2moPoV^^ and m* is the ef- 
fective mass given by m* = too(7c + '^Ep / (ZUih) + 
Ep/{Wso))'\ Uih = E-Hih, Uso = E-Hso, and Hih 
and Hso are the diagonal elements of the light-hole and 
spin-orbit split-off bands in the eight-band Hamiltonian. 
From Eq. we find that the dominant contribution 

to the RSS consists of the interface term AE'i^'' in the 
valence band and the external electric field term /S.En ' , 



Ai?,«(fc||) = — fc|| ^ |F„(z,)l' [7(0 - 7(^7)]' (2) 
i 
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EpeFk\\ I dz \F„iz)\'{U,l'^U,-'). 

(3) 
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Here Fn{z) is the envelope function of the nth subband 

along the z axis, and {zj} denote the z coordinates of 

the interfaces. According to the Ehrenfest's theorem, 

(2) 

the external electric field contribution AEn is approxi- 
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FIG. 1: (Color online) Comparison of the linear (red or 
dark gray lines) and nonlinear (black lines) models with 
the numerical results (green or light gray lines) obtained 
from the eight-band k ■ p theory for the lowest conduc- 
tion subband of 15-nm-wide (a) GaAs/Gao.eTAIo.ssAs, (b) 
Gao.47Ino.53As/Alo.48lno.52As, (c) InAs/Ino.75Alo.25As, and 
(d) Hgo.79Cdo.21Te/CdTe QW's with an electric field 60 
kV/cm. 



mately cancelled by the interface electric field contribu- 
tion in the conduction band.[28j As a result, the domi- 
nant contribution to RSS comes from the interface term 
only.*^ From Eq. ||2Jl, we see that RSS is approxi- 
mately inversely proportional to the effective band gap 
Eg-f^ = E ~ Ey, where E^ is the valence band edge 
and E is the electron energy. Since the electron energy 
E increases approximately quadratically as fc|| increases, 
the RSS would always begin to decrease after has ex- 
ceeded a critical value. This reveals that the previously 
found deviation of the RSS from linear behavior in cer- 
tain QW'sUl m m is actually a universal behavior, 
suggesting that the widely used linear Rashba model need 
to be checked and, if necessary, replaced by a nonlinear 
one. 

Enlightened by above discussions, we propose the fol- 
lowing two- CO efficient nonlinear Rashba model to de- 
scribe the RSS of a given subband n. 



A£;„(fc||) 



2a„fc|| 
1 + I3nkr 



(4) 



where q;„ is the linear Rashba coefficient of the nth sub- 
band for the widely used linear Rashba model /S.En (fc|| ) = 
2a„fc||, while /3„fc^ describes the contribution from the 
kinetic energy of electron in the n-th subband. The 
latter leads to the decrease of RSS when fcy exceeds a 
critical value feg = l/^fiJ^, which is determined from 

['^(Ai^"(fc||))/*ll]fc,=,o=0- 

In Fig. 1 , we compare the linear and nonlinear Rashba 
models with the numerical results obtained from the 
eight-band k • p theory for different QW's. The band 



FIG. 2: Rashba coefficient a/F (upper panels) and 
P (down panels) of the lowest conduction subband 
as a function of the well width and external electric 
field for Gao.47Ino.53As/Alo.48lno.52As [(a) and (b)], and 
Hgo.79Cdo.21Te/CdTe [(c) and (d)] QW's. The units of a/F 
and P are meV nm/(100 kV/cm) and nm^, respectively. 



gap of the well material ranges from 1.519 eV (for GaAs) 
to 0.1 eV (for Hgo.79Cdo.21Te). We see that the deviation 
of the linear model from the numerical results could be 
surprisingly large, especially for narrow gap QW's, while 
the nonlinear model is in good agreement with the nu- 
merical results. Thus we substantiate the necessity to use 
the nonlinear Rashba model instead of the linear one. 

To characterize the nonlinear Rashba model, we con- 
sider the dependence of the Rashba coefficients a„ and /3„ 
on the various band parameters. To the lowest order, the 
probability asymmetry at the two interfaces of the QW 
is proportional to the external electric field F, then Eq. 
@ suggests AS„(k||) ^ k\\F/[Eg + E^o + h^k^^/{2m;j], 
where Eno is the quantum confining energy and to* is 
the effective mass of the nth subband. This shows that 



F 



En + EnO ' 



f3n 



m*{Eg + Eno) ' 



(5) 
(6) 



a„ and /3„ both increase with decreasing band gap and 
subband index n, or increasing well width, i.e, the quan- 
tum confining energy Eno- The significant difference is 
that an increases with increasing external electric field 
F, while Pn is approximately independent of F. No- 
tice, however, additional dependence of a„ and /?„ on 
the external electric field may comes from the quantum 
confining energy, especially at large electric fields. 

In Fig. 2, we plot the Rashba coefficient a/F and /3 
of the lowest conduction subband (obtained by fitting 
the results from the eight-band k • p theory) as functions 
of the well width and external electric field for differ- 
ent semiconductor QW's. Some common features can be 
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observed. Firstly, a/F and (3 are weakly dependent on 
the external electric field. This shows that a is approx- 
imately proportional to the external electric field, while 
(3 is approximately independent of F, in agreement with 
Eqs. Q and lO- Secondly, a/F and (3 increase with 
the increase of the well width and begin to saturate at 
large well width. The saturation value decreases with 
the increase of the electric field, and this trend becomes 
increasingly pronounced when the band gap of the well 
decreases. The increase and saturation behavior comes 
from the competition between the band gap Eg and the 
quantum confining energy EnQ especially at small well 
width. The decrease of the saturation values of a/F and 
(3 at large well width comes from the enhancement of 
the quantum confining energy Eno due to the triangular 
potential induced by the external electric field. It also re- 
flects the fact that the relationship a (x F begins to break 
down for wide QW's. Thirdly, the critical well width at 
which a/F begins to saturate decreases with the increase 
of the electric field. This can be understood from the 
competition between the square potential [produced by 
the conduction band profile Ec{z)\ and the triangular po- 
tential (produced by the external electric field) . With the 
increase of the well width, the square potential becomes 
weaker and the triangular potential begins to dominate. 
This transition occurs at a smaller well width when the 
external electric field get stronger, leading to the satu- 
ration of the quantum confining energy Eno and, conse- 
quently, the Rashba coefficients a/F and /? at a smaller 
well width. The similarities between a/F and /? are in 
agreement with our analytical discussions in Eqs. (jSJ and 
However, we should notice that a is approximately 
linearly dependent on the external electric field F , while 
(3 is approximately independent of F . 

In the above, we see that the nonlinear Rashba model 
differs significantly from the linear model. As a result, 
we expect that this nonlinear Rashba model would lead 
to a series of modifications to the various spin-related 
properties of the electron, which have been inve stig ated 
based on the linear Rashba model. [13, El, EE EE 113, 
EEEEEEimilllEil Asan example, we consider 
the electron spin relaxation caused by Rashba spin-orbit 
coupling in two-dimensional electron gas. The electron 
spin relaxation process has received intensive interest re- 
cently because it plays an essential role in the practical 
application of s pint ronic devices and quantum informa- 
tion processing. |3ll [sE HE 1^ For electron in n-doped 
semiconductors at large electron density, DP mechanism 
is the dominant pro cess and it has been investigated 
by several groupsjisl ll^ using the linear Rashba 
model in the framework of the density matrix formal- 
ism. We follow this procedure and calculate the DP 
SRR [F,, = Tyy = r,,/2, F„- = (z j)] using the 
eight-band k • p theory, the nonlinear and linear Rashba 
models, respectively. The results for different QW's are 
shown in Fig. 3. The SRR's from the nonlinear Rashba 
model show pronounced suppression at large Fermi wave 
vector kp (i.e., electron density), in good agreement with 
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FIG. 3: (Color online) Low temperature DP SRR in the 
lowest conduction subband of 20nm- and lOnm-wide (a) 
GaAs/Gao.erAlo.ssAs, (b) Gao.47Ino.53As/Alo.48lno.52As, (c) 
InAs/Ino.75Alo.25As, and (d) Hgo.79Cdo.21Te/CdTe QW's 
with an electric field 60 kV/cm and momentum relaxation 
time 0.1 ps. The results are obtained from the eight-band 
model (green or light gray lines), the nonlinear (black lines) 
and linear (red or dark gray lines) Rashba models, respec- 
tively. 
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FIG. 4: Tzz (in units of ps) for 10-nm-wide 
InAs/Ino.75Alo.25As [(a) and (b)] and Hgo.79Cdo.21Te/CdTe 
[(c) and (d)] QW's with a typical momentum relaxation time 
0.1 ps. The linear (nonlinear) Rashba model is used for (a) 
and (c) [(b) and (d)]. 



the numerical results from the eight-band k • p theory, 
while those from the linear Rashba model increase mono- 
tonically with fc^?, which qualitatively disagrees with the 
nonlinear model. The difference between the linear and 
nonlinear Rashba models increases significantly with the 
increase of fci?. Taking the 20 nm wide QW's for exam- 
ple, at kp — 0.4 nm^^ (corresponding electron density 
Uf, K, 2.5 X 10^^ cm~^, which is of interest to, e.g., trans- 
port measurements), the relative difference between the 
two models reaches 40%, 140%, 500%, and 3800% for 
Figs. 3(a)-3(d), respectively. 
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In the above, we consider only the Rasbha spin-orbit 
couphng in order to demonstrate the significant conse- 
quences that are introduced by our nonhnear Rasbha 
model. To perform realistic calculations, however, both 
the nonlinear Rashba model and the Dresselhaus spin- 
orbit coupling should be included. Recently it was pre- 
dicted based on the linear Rashba model that the in- 
terplay between the Rashba and Dresselhaus spin-orbit 
interactions would greatly enhance the electron spin life- 
time in [lll]-orientated quantum wells when the linear 
Rashba coefficient a and the linear Dresselhaus coeffi- 
cient fin are opposite (i.e., when the ratio ~a//3D = 
1)-E3 As a second example, we reconsider this problem 
using the nonlinear Rashba model. The electron spin life- 
time Tzz (along the growth direction of the QW) obtained 
from the two models are compared in Fig. 4 as functions 
of the ratio (—a/f3]j) and Fermi wave vector kp. The 
black regions correspond to the resonant enhancement 
regions of the spin lifetime. The most striking feature 
is that the ratio (— a/zSu) at the resonance peak (latter 
referred to as the peak position) shows non-motonic de- 
pendence on the Fermi wave vector [cf. Figs. 4(b) and 
4(d)], which differs qualitatively from the monotonically 
decreasing behavior of the linear Rashba model [cf. Figs. 
4(a) and 4(c)]. Here the monotonic decrease of the peak 
position with increasing kp in the linear Rashba model 
comes from the Dresselhaus k^ term, while the nonlin- 



earity of the Rashba effect competes against the Dressel- 
haus k^ term, leading to non-monotonic dependence of 
the peak position on the Fermi wave vector. The values 
of the spin lifetime obtained from the nonlinear Rashba 
model also differ from the prediction of the linear Rashba 
model by up to several orders of magnitude [cf. Figs. 4(c) 
and 4(d)]. 

In summary, we have revealed that the RSS in semicon- 
ductor QW's is intrinsically a nonlinear function of the 
wave vector. It may be overestimated significantly by the 
linear Rashba model, especially in narrow-gap QW's. We 
propose a two-coefficient nonlinear Rashba model, which 
is in good agreement with the numerical results obtained 
from the eight-band k • p theory. Using this nonlinear 
model, we found pronounced suppression of the DP SRR 
at large electron density, and a non-monotonic depen- 
dence of the ratio {—a/(3r)) at the resonance peak of the 
electron spin lifetime on the Fermi wave vector in [111]- 
oriented QW's, both in qualitative disagreement with the 
predictions of the linear Rashba model. The values of the 
spin lifetime obtained from the two models may differ by 
up to several orders of magnitude. Further surprising re- 
sults are expected when the nonlinear Rashba model is 
applied to other fields. 
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